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1>.  AStTAACT 

The  turbulent  diffusion  of  a  pa 

ssive  scalar  in  a  stratified  fluid 

is  studied  using  a  simplified  model  for  the  calculations.  A  semi-empirical 
formula  for  the  diffusion  coefficient  is  obtained  as  a  function  of  turbulent 
energy  dissipation  rate  and  scale  of  the  diffusing  patch.  The  formula  is 
demonstrated  to  give  values  which  are  consistent  with  known  experimental 
results.  It  is  suggested  that  the  diffusion  coefficient  at  intermediate  depths 
is  an  order  of  magnitude  smaller  than  that  observed  in  near  surface  waters. 
Using  these  values  of  the  diffusion  coefficients,  calculations  show  that  the 
mean  concentration  for  a  patch  initially  100 meters  across  decreases  by  a 
factor  of  10  in  two  hours  near  the  surface,  but  in  eight  hours  at  intermediate 
depths. 
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Abstract 


The  turbulent  diffusion  of  a  passive  scalar  in  a  stratified  fluid  is  studied 
using  a  simplified  model  for  the  calculations.  A  semi-empirical  formula  for 
the  diffusion  coefficient  is  obtained  as  a  function  of  turbulent  energy  dissipation 
rate  and  scale  of  the  diffusing  patch.  The  formula  is  demonstrated  to  give 
values  which  are  consistent  with  kr.^wn  experimental  results.  It  Is  suggested 
that  the  diffusion  coefficient  at  intermediate  depths  is  an  order  of  magnitude 
smaller  than  that  observed  in  near  surface  waters.  Using  these  values  of  the 
diffusion  coefficients,  calculations  show  that  the  mean  concentration  for  a  patch 
initially  100  meters  across  decreases  by  a  factor  of  10  in  2  hours  near  the  sur¬ 
face,  but  in  8  hours  at  intermediate  depths. 
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1.0  Introduction 


The  turbulent  diffusion  of  a  passive  scalar  in  a  stratified  fluid  was 
studied  using  three  models  for  turbulent  difiusivity.  The  calculations  are 
approximate  in  the  sense  that  only  diffusion  in  one  horizontal  direction  is  con¬ 
sidered  and  that  the  effects  of  the  turbulence  are  "modeled".  However,  the 
"model"  diffusion  equations  are  solved  exactly,  in  some  cases,  or  numerically 
with  high  accuracy  in  all  other  cases. 


It  is  assumed  that  the  quantity  whose  diffusion  is  being  studied  is  a 
passive  scalar.  That  is,  the  quantity  undergoing  diffusion  can  be,  like  temp¬ 
erature  or  concentration  of  a  substance,  specified  by  a  scalar  function  of  space 
and  time  and  the  dynamics  of  the  flow  are  totally  unaffected  by  the  variation 
in  concentration  of  the  diffusing  quantity.  It  is  further  assumed  that,  what¬ 
ever  the  diffusing  substance  is ,  it  does  not  undergo  any  chemical  reaction 
with  significant  release  or  absorption  of  heat. 

Finally  it  will  be  assumed  that  the  fluid  is  stratified.  In  many  areas 
of  the  ocean  there  is  an  appreciable  variation  of  density  with  depth  outside  of 
a  well  mixed  surface  layer.  It  has  been  known  for  many  years  (see  Sverdrup,  et  al 
1942  )  that  the  presence  of  a  stable  stratification  reduces  the  turbulent  inten¬ 
sity  and  inhibits  vertical  transport  of  both  scalar  and  vector  properties „ 

Okubo  and  Pritchard  (1969)  give  measured  values  of  the  horizontal  and  vertical 


diffusion  coefficients  lor  the  ocean.  The  horizcntal  diffusion  coefficients  are 


3  3?  + 

in  the  range  of  10  to  5  x  10  cm  /sec  whi»e  the  vertical  diffusion  coefficients 


+  Ozmidov  (1969)  quotes  values  ranging  from  10**  to  10^cm^/sec 
for  the  coefficient  of  horizontal  diffusivity. 
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in  the  upper  mixed  layer 
in  a  significa  it  thermocline 
in  the  deep  layer 
near  bottom 


10  - 100  cm^/sec 
0.  01  -  1.  Octn^/sec 
0.  1  -10.  0  cm2/  sec 
1-10  cm^/sec. 


It  is  clear  that,  excluding  the  upper  mixed  layer,  whichever  of  these  values 
is  taken  for  the  inefficient  of  vertical  diffusivity,  the  ratio  of  the  coefficient 
of  horizontal  diffusivity  to  the  coefficient  of  vertical  diffusivity  is  greater  than 
100.  It  therefore  seems  r^ason^He  to  neglect  the  vertical  diff  ision  and  to 
consider  only  diffusion  in  the  horizontal  plane. 

The  source  of  tie  scalar  is  assumed  to  be  a  self-propelled  body 
traveling  through  the  ocean.  The  scalar  is  assumed  to  be  released  into  the 
wake.  The  wake  grows  by  entrainment  and  eventually  collapses  (see  the  model 
of  Ko,  1971).  If  the  scalar  is  uniformly  distributed  throughout  the  wake  just 
behind  the  body,  the  model  of  the  growth  and  collapse  of  the  wake  describes  the 
spread  of  the  scalar.  That  is,  the  scalar  is  uniformly  (on  average)  distributed 
within  the  wake  and  the  calculation  of  the  wake  boundary  as  a  function  of  time 
gives  the  mean  distribution  of  the  scalar  as  a  function  of  time.  The  more  dif¬ 
ficult  question  of  the  diffusion  of  the  scalar  from,  say,  a  point  source  within 
the  wake  will  be  considered  in  a  following  report. 

Ko's  model  shows  that  the  turbulent  intensity  decreases  rapidly  with  time 
after  the  beginning  of  collpase.  When  the  turbulent  intensity  within  the  wake 
drops  appreciably  below  the  ambient  turbulent  intensity  of  the  ocean,  Ko's  model 
can  no  longer  describe  the  diffusion  process  and  the  u-rbulence  cf  the  ocean 
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controls  the  diffusion  process. 

In  the  calculation  reported  here,  the  "equivalence  principle"  is  used, 
i.  e. ,  the  model  is  that  of  a  "slice  of  the  ocean"  in  a  vertical  plane  perpendicular 
to  the  axis  of  the  wake.  Since,  as  discussed  above,  vertical  diffusion  can  be 
neglected  compared  to  horizontal  diffusion,  the  problem  can  be  reduced  in  this 
approximation  to  that  cf  turbulent  diffusion  in  one  space  dimension  and  time. 

In  fact,  the  problem  is  really  a  function  of  both  spatial  coordinates  in 
the  horizontal  plane  and  of  time.  Within  the  approximation  used  in  this  report, 
for  example,  the  meander  of  the  wake  due  to  large  eddies  cannot  be  described. 
This  is  a  serious  limitation.  A  possible  experiment  would  involve  sampling 
the  wake  along  some  path.  !f  the  sample  were  taken  along  a  path  parallel  to 
the  track  of  the  body,  the  meanders  of  the  wake  could  cause  large  variations  in 
the  observed  concentration.  The  calculations  reported  here  are  of  the  average, 
over  an  ensemble,  of  the  concentration.  While  the  variation  of  the  mean  con¬ 
centration  as  a  function  of  time  is  important,  it  is  not  the  whole  story.  The 
variation  of  the  concentration  with  time  in  horizontal  plane  is  now  being 
investigated  and  will  be  reported  at  a  later  date. 


f 
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2.  0  Formulation  of  the  Problem 


Let  be-  the  spatial  coordinate  parallel  to  the  fi-ee  surface  and 
perpendicular  to  the  path  of  the  self-propelled  boiy,  be  the  time  and 
S  be  the  mean  concentration  of  some  passive  scalar.  The  concentration 
S ( then  governed  by  a  diffusion  equation 


dS  __  e) 

H-  ST) T(  A  dK 


*s). 


(') 


It  has  been  assumed  that  mean  velocity  of  the  fluid  is  zero.  ^  is  the 
effective  diffusion  coefficient.  If  the  fluid  is  *rat  rest1',  i.e.,  there  is  no 
turbulence,  Jie.n  is  the  molecular  diffusivity,  but  if,  as  is  usually  the 
case,  there  is  turbulence,  theu^^  is  an  eddy  diffusivity. + 

It  is  assumed  that  the  concentration  is  uniform  in  a  region  of  width  y 
at  some  initial  time-  which  is  taken  to  be  ^  -  0.  As  discussed  in  the  intro 

duction,  this  is  some  arbitrary  time  after  the  wake  has  collapsed  and  the  in 
tensity  of  the  wake  turbulence  has  dropped  below  the  ambient  turbulence  in¬ 
tensity.  Thus  th2  initial  condition  is  that 


S~St  j  f  =  o  ,  -J./z  ty&J'/i  .  (i) 


+ 


See  Appendix  A  for  a  detailed  discussion  of  the  derivation  of  equation  (1) 
and  the  approximations  used. 
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Without  any  loss  of  generality  S*#  can  he  taken  equal  to  1 , 
Equation  (1)  can  be  put  into  dimensionless  form  by  defining  appropriate 
length  and  time  scales.  Let  the  length  scale  be  A  and  let  K  be  a 
reference  diffusivity  (dimension  (length)^ /time).  The  time  scale,  ,  is 

then 


A  »  A  VX  .  f» 


«r* 

1 


I 


Define  the  dimensionless  spatial  coordinate 


f  =  , 


(vj 


and  the  dimensionless  time  coordinate 


J  -  */*.  • 


(sJ 


Under  this  transformation,  equation  (I)  becomes 


J_f  JS  ) 

d'h  ~  df  [A*  If  J 

where  the  dimensionless  diffusion  coefficient 


TfJ  «■> 


S"  -  A' A".  . 


(7  J 


The  initial  condition  is  now, 

S=1  ;  7*0,  -i/x  'A-  (V 

The  problem  is  uniquely  specified  whei^^/^f  { ^  is  given.  Since 
a  general  olution  to  the  turbulence  problem  is  not  available  ,yti  must  be 
approximated  or,  better,  modeled.  A  common  approximation  is  to  use  a 
constant  value  for  the  eddy  diffusivity.  This  is,  at  best,  a  very  rough  approxi¬ 
mation  and,  of  course,  is  crucially  dependent  on  the  choice  of  a  numerical 
value  of •  The  value  sf^^/  not  only  varies  from  one  location  to  another 
in  the  ocean  but  also  varies,  at  a  given  location,  with  the  scale  of  the  diffusing 
patch.  Because  the  constant  eddy  diffusivity  model  is  so  widely  used,  and  be¬ 
cause  other,  more  appropriate  diffusion  models  can  sometimes  be  reduced, 
under  a  coordinate  transformation,  to  an  equation  with  a  constant  eddy  diffusi¬ 
vity,  this  model  will  be  treated  first. 

The  variation  of (or  *  )  with  patch  size  can  be  understood  in 
terms  of  the  variation  of  energy  in  the  spectium  of.  turbulent  eddies.  Consider 
a  patch  of  size  .  Turbulent  eddies  which  are  much  larger  than  simply 
convect  the  patch  without  appreciably  distorting  cr  spending  it.  Turbulent 
eddies  which  are  much  smaller  than  ^  stir  up  the  patch  and  provide  the  fine 
scale  mixing.  Of  course  the  small  eddies  do  prrvice  enhanced  mixing  at  the 
edges  of  the  patch.  However,  it  is  the  eddies  v  ith  a  seals  approximately  equal 
to  that  provide  the  major  mixing.  These  eddies  distort  the  patch  and 


greatly  increase  the  surface  area  of  the  patch  and  thus  enhance  the  mixing  due 


to  the  small  eddies  at  the  edges.  Clearly  the  most  important  quantity  is 

the  rate  of  distortion  of  patch  by  the  eddies  ox  scale  J  .  But,  the  turbulent 

energy  Mit X  increases  with  decreasing  1*1  ,  i.o.,  with  increasing 

scale,  at  least  up  to  the  maximum  scale  of  the  turbulence.  Therefore,  the 

larger  the  patch,  the  more  energy  there  is  in  the  turbulent  eddies  of  the 

same  scale  which  distort  the  patch.  Of  course  there  is  an  upper  limit  in 

that  the  turbulent  energy  spectrum  does  not  increase  indefinitely  with  decreasing 


III  •• 


there  is  a  peak  in  the  spectrum  at  the  scales  at  which  the  turbulence 


is  Leing  generated.  This  is,  however,  unimportant  in  this  problem.  The  peak 
in  the  spectrum  lies  at  scales  greater  than  tens  of  kilometers  and  the  patches 
to  be  studied  here  have  sizes  in  range  1  meter  to  10  kilometers. 

Richardson  (1 926)  studied  the  turbulent  dispersion  of  particles  in  the 
atmosphere  and  showed  experimentally  that 


k  = 


(V 


^  2/3  §  f 

where  £  is  a  constant  (dimensions  cm  /sec  when  /{,  is  in  cm)  and  sC 

is  the  separation  of  the  particles.  A  number  ox  other  experiments  have  been 
performed  on  v/ater  surfaces,  Richardson  &  Stommel  (1948),  Stommel  & 
Marin  (1949),  Ozmidov  (1959)  and  Gray  &  Pochapsky  (1964),  and  have  also 
shown  agreement  with  equation  (9).  In  a  forthcoming  report,  Kuo  (1972)  will 
discuss  a  theoretical  justification  of  the  "4/3"  law.  In  view  of  this  experi¬ 
mental  evidence  and  theoretical  justification,  it  will  be  assumed  that  the  best 
model  of  the  turbulent  diffusion  coefficient  is  the  "4/3"  law  (equation  (9)). 
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The  remaining  conceptual  problems  are  (1)  what  is  the  appropriate  scale, 
i.e.,  what  is  ^  and  (2)  what  is  the  value  of  £  ?  These  questions  will 

be  discussed  in  detail  below.  First,  however,  a  number  of  eddy  diffusivj.ty 


3.  0  Diffusivity  Models 


In  this  section,  three  models  for  the  eddy  diffusivity  are  examined, 
namely,  a  constant,  a  variable,  and  an  integral  scale  value.  The  resulting 
mean  concentration  is  calculated  for  each  of  these  models  in  terms  of  a 
dimensionless  time  variable. 
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I 

I 

I 

I 

I 


n 


3. 1  Constant  Eddy  Diffusivity  Model 


If  the  eddy  diffusivity  is  constant. 


Equation  (6)  is  then 


c>S  _  ys 

dy  ~  df ' 

with  the  initial  condition  given  hy  (8) . 
The  solution  is  easily  found  to  be 

with  er-f(i)  defined,  as  usual, 

"Z  f 

=  -p=,  j  e 
nr  j 
o 


„  U  =  (ro) 


j 


ov 


The  three  quantities  of  interest  are  (1)  the  spatial  distribution  of 
the  scalar  at  a  fixed  time  (given  by  equation  (12));  (2)  the  concentration  at 
the  center  of  the  patch  as  a  function  of  time 


1 
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1 


1 

’6 


:| 

£5 


"S 
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& 


nrtig^l#  'Z'-'^-rTC- -I's = * 


**  /r-  ^ 

Then,  if  «^C  is  proportional  to  ^  , 


A /A  - /OV77  =  fi+Wp t  ( if) 

where  the  zero  subscript  refers  to  values  at  *^  =  0  (  ^  =0). 

Equation  (15),  using  equations  (12)  and  (14),  has  been  solved  numeri¬ 
cally  for  Q?  =0.1  and  0.25.  The  results,  curves  of  A /A  as  a  function 
of  ir/tt  are  plotted  in  Figure  (1),  along  with  values  of  the  scale  defined  in 
terms  of  the  variance,  i.e.,  A/A  from  equation  (18). 

It  can  be  seen  from  the  results  plotted  in  Figure  (1)  thr-t  there  is  little 
difference  between  the  size  of  the  patch  as  given  by  the  different  definitions  of 
AfA  .  As  expected,  the  largest  size  is  obtained  by  defining  J  as  the 
1/10  concentration  point.  The  scale  ..izes  defined  by  the  1/4  concentration 
point  and  the  variance  agree  within  about  10%  and  are  about  25%  smaller  than 
the  scale  defined  by  the  1/10  concentration  point.  These  scales  are  the  same 
to  within  the  accuracy  of  the  approximations  which  are  imbedded  in  the  calcu¬ 
lations.  In  particular,  there  is  uncertainty  in  the  choice  of  an  eddy  diffusion 
model  and  error  in  the  use  of  an  eddy  diffusion  to  represent  the  turbulence  in 
the  first  instance.  Finally,  the  scale,  as  calculated  here  is  an  ensemble 
average  and  in  any  single  realization  the  actual  scale  will  differ  from  the  en¬ 
semble  average.  In  view  of  these  uncertainties,  it  seems  reasonable  to  use 
the  variance  definition  of  the  scale  as  a  reasonable  approximation.  The  mean 
concentration  is  then  just  the  reciprocal  of  A/A.  .... 

<S>  *  [/+i*f  +/U)]~  .  (it) 
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3.2  Variable  Eddy  Diffusivity  Model 


One  possible  choice  for  the  scale  to  be  used  in  the  "4/3"  law  is  to 


take 


This  model  is  based  on  the  idea  that  the  scalar  is  diffusing  away  from  the 
center  of  the  patch  and  therefore  the  appropriate  scale  is  the  distance  from 
the  center  of  the  patch.  With  this  approximation 


-  2  V 


ft  0 


and 


X  =  ?  A ih. 


(it) 


Therefore 


t*3J 


and  equation  (6)  becomes 


c)  I  c>S 
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It  can  be  shown  using  elementary  techniques,  see  Appendix  B,  that 
the  solution  to  equation  (24)  which  satisfies  the  initial  condition  (8),  is 

S(?s  J)  ~  "t  j  et~f( zr)  - 


2 

nr 


[ - i. e-B'l J  (is) 


*i)*/f^  .  (17) 


The  concentration  at  the  center  of  the  patch  is , 


5(°,  •Jj  -  e,'f'(B)~'f=pie  *,  (20 

i  ■=  V/2  "h  Tf  'h)  _ 


with 


iiy'nmMfywprit  i*1 


If  the  scale,  or  size,  of  the  patch  is  determined  from 

sUAy***^?),  a*) 

then  equaiicnc  (25)  through  (30)  can  be  solved  numerically  for  J/A 
function  of  ^  ^ ^/^C^OT  an'/  particular  choice  of  CL*  •  If,  however 


[for  any  particular  choice  of  CL  •  If.  however,  the 


scale  is  defined,  as  before,  in  terms  of  the  variance,  it  can  be  shown  (Appen¬ 


dix  B)  that 


rh>  °TT  +  (Wr)?, 


no 


and  thus  that 


-*H  * h+V8rh’‘ .  („) 


The  scale  size 


j/j. 


,  as  determined  from  equation  (30),  is  plotted 


in  Figure  (2)  as  a  function  of 


i/i. 


for  a  =0.1  and  0.25  as  well  as  the 


scale  size  defined  in  terms  of  the  variance,  equation  (32).  From  figure  (2) 
it  is  seen  that  the  scale  is  approximately  equal  to  until  C  ZZ  0.02 

and  then  begins  to  grow  rapidly. 

In. contrast  to  the  constant  eddy  diffusivity  model,  the  variance  scale 
is  'ei'ger  than  0.10  and  0.25  concentration  scales.  This  can  be  understood  in 
terms  of  the  eddy  diffusivity  model,  equation  (23).  The  rate  of  diffusion  in¬ 
creases  with  increasing  distance  from  the  center  of  the  patch.  This  tends  to 
give  slow  diffusion  near  the  center  of  the  patch  and  fast  diffusion  near  the  edges. 
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This  implies  that  there  is  a  much  sharper  "edge"  to  the  patch  than  in  the 
ca.'e  of  a  constant  eddy  diffusivity  and  that  there  will  be  a  long  "tongue"  of 
material  at  low  concentration  extending  out  from  tie  edge.  Since  the  material 
in  the  "tongue"  is  diffusing  rapidly  compared  to  the  material  at  the  center,  the 
"tongue"  grows  rapidly.  Now  in  calculating  the  variance  the  material  at  large 
^  is  heavily  weighted  (by  j?  ^  )  and  thus  contributes  more  than  material 
near  the  center.  Because  of  the  sharp  "edge"  and  long  "tongue",  a  large  dif¬ 
ference  between  variance  scale  and  constant  concentration  scale  is  to  be  ex¬ 
pected  . 

A  reasonable  approximation  in  this  case  would  appear  to  be  to  take 
the  scale  as  the  variance  scale.  The  .mean  concentration  is  then 


-'A 


<s>  •[«■(?£)(*/+?]  .& 
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3.3  Integral  Scale  Eddy  Diffnsivity  Model 


Finally,  it  can  be  argued  that  the  appropriate  scale  for  the  diffusion 
is  constant  at  any  instant  of  time  and  is,  in  fact,  the  "size"  of  the  patch. 
Taking 


and 


X  * 


t 

(3VJ 

j* 

j 

C3S-J 

t/M  )  4/3 1 

f, 

(3<>) 

=  **5(*A)t 

(3V 

or  perhaps  in  terms  of  the  variance  scale, 


-V-4  =  fo~Vrp '  .  f 


The  mathematical  problem  is  now  to  solve 

<=>S  ,  t  <?*5 

d*h  ~^(V 


J 


fl  ,  (39j 


S (f,o)*  1  j 


-19- 


:»iaatatatias»a8>»fe>a?»««a?e<>^»i 


. ■  HIPJI  i  tun  n  'I1  T  mnFrwitFirwinirin--rT — nrm-—— — a—,v 


Since^t/  is  only  a  function  of  ,  it  is  possible  to  eliminate 
from  equation  (39)  by  a  suitable  coordinate  transformation.  De¬ 
fine  the  weighted  time  coordinate  by 


M 


c>  d7~  c)  .  t  £ 

<*J  c>7{  e>T  ^  ^  JTj 


and  equation  (39)  becomes 


i£ 

c>T" 


£s_ 

or 


with  the  same  initial  condition. 


(V*) 


0/3) 


This  is,  of  course,  the  diffusion  equation  for  a  constant  eddy  diffusivity. 
The  solution  is  equation  (12)  with  ^  replaced  by  '£“■  .  In  order  to  interpret 
the  results  it  is  necessary  to  transform  back  from  the  to  the  0^  variable. 

(Recall  that  ^  is  directly  proportional  to  the  time  while  is  a  weighted 

integral  of  the  time.)  The  diffusion  scale  (  Jf/A  )  can  be  calculated,  whatever 
its  definition,  as  a  function  of  .  Then  the  eddy  diffusivity  can  be  calculated 

from  (  *14  ).  NoWjthe  edd>  diffusivity  must  be  the  same  in  the  0^  and  ^ 
coordinates  when  it  is  calculated  at  corresponding  (not  equal)  values  of  ^ 


.  The  correspondence  is  that  determined  by  the  transformation. 
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equation  (41).  From  the  differential  form  of  this  transformation 


cl.'i-  -  yt  (y)  cl  %  t 


(yy) 


can  be  obtained  by  using  the  above  principle.  Thus 


sytd,)  , 


Ns) 


and  correspondine  %  and 


,  so  that 


c(t 


u  c  / 

-  =  cl 

yi't)  ( 


(Yl) 
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~  J  yffTr') 


(yV 


Using  a  constant  concentration  to  define  the  scale,  values  of  A/A 
as  a  function  of  2"  are  known  from  the  numerical  solution  of  equations  (12), 
(14)  and  (15).  The:  ing  equations  (36)  and  (4?)  a  numerical  integration 

yields  ^(z)  .  This  calculation  has  been  carried  out  with  Oj*  (the  value  of 
the  concentration  defining  the  edge  of  the  patch)  equal  to  0.1  and  0.25.  Since 
(A /A  )  and  S  ( o,  T)  are  known  as  functions  of  ,  it  is  easy  to  tabu¬ 

late  (  A/A  > as  *  £u"clion  ot  ^  and  also  S  (  cj  ^)- 

On  the  other  hand,  if  the  scale  is  defined  in  terms  of  the  variance,  the 
transformation  can  be  inverted  analytically.  In  terms  of  the  variable. 


* 

i 


Figure  (3)  is  a  plot  of  {  Jl/J0  )  as  a  function  of  (  -&/-60  )  for  the 
constant  concentration  scales,  £^'*’ =  0.1  and  0.25,  and  the  variance  scale. 
Again  the  largest  scale  is  that  obtained  from  the  1/10  concentration  point. 
The  variance  scale  and  the  1/4  concentration  point  scale  are  very  close. 


within  10%.  It  appears  to  be  reasonable  to  use  the  variance  definition  of  the 
scale.  The  mean  concentration  is  the  leciprocal  of  JMo  and  is 

<S>  =  /  /+*(*/+)]  (S3) 


css  patch  si'/.e  (  ./  /j\  4  as  a  function  of  the  dimensionless  time 
1  4/3  law  (  r-.  |  (i  /  ^  J  */'■* )  has  been  used.  The  different  curves 

definitions  of  the  patch  size  discussed  in  the  text. 


4.0  Interpretation  of  the  Numerical  Results 


4.1  Theoretical  Cons idt  rations 


Although  the  constant  eddy  diffusivity  model  has  often  been  used,  it 
is  clear,  on  both  experimenal  and  theoretical  grounds,  that  the  "4/3  law" 
provides  a  better  model  of  the  diffusion  process.  While  choosing  a  scale 
proportional  to  the  distance  from  the  center  of  the  patch  (a  local  scale)  may 
be  a  reasonable  approximation  for  small  times,  the  integral  scale  is  prob¬ 
ably  the  most  reasonable  scale  to  use  in  the  4/3  law.  Figure  (4)  shows  the 
variation  of  concentration  at  the  center  of  the  patch  %  J  and  the  mtar. 

concentration  <*>  with  time  Ufa  )  for  all  three  models.  From  now  on 
all  the  discussion  will  center  on  the  integral  scale  4/3  law  Mf 

In  order  to  interpret  these  results,  the  initial  length  scale,  , 
and  initial  time  scale  to  must  be  determined.  There  is  no  problem  with 
the  length  scale;  should  be  the  initial  size  of  the  diffusing  patch.  As 


to 


t.  *  -f.'/x.  , 


fry) 


the  initial  eddy  diffusivity 


X  -  rJ. f/’  , 


( ss) 


is  required  for  definition.  In  order  to  pick  a  value  of  £  ,  experimental 

data  must  be  used. 
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Figure  4.  Variation  of  the  concentration  at  the  center  of  the  patch 
(dashed  line)  and  the  mean  concentration  (solid  line)  as 
a  function  of  the  dimensionless  time  (  t/i9).  The  curves 
as  labeled  correspond  to  the  three  different  diffusion  laws 
discussed  in  the  text. 


Most  experimental  studies  of  horizontal  diffusion  in  the  ocean  have 
been  conducted  in  near  surface  waters,  if  not  on  the  surface,  see  Ozmidov 


(1959)  and  Okubo  &  Pritchard  (1969).  Although  there  is  considerable  scatter 
in  the  data,  a  reasonable  fit  appears  to  be  given  by  equation  (55)  with 

^ 5*  0.01  cm^^/sec  (  jP  in  cm). 

For  the  present  purposes  these  measurements  near  the  surface  are 
not  satisfactory.  The  near  surface  waters  are  disturbed  by  wind  driven 
currents  and  surface  waves,  for  example,  and  generally  are  not  stratified. 

It  would  be  desirable  to  have  an  extensive  set  of  measurements  of  the  diffu- 

0  00 

sion  coeffcients  as  a  function  of  scale  for  various  depths  and  Vaisala  periods. 
Such  experimental  results  are  not  available.  There  appears  to  be  only  a  few 
diffusion  measurements  at  depth.  Fortunately,  some  additional  measure¬ 
ments  are  available,  which,  properly  interpreted,  will  give  further  informa¬ 
tion  on  the  variation  of  ^  with  depth  (  0  )  and  Vaisala  frequency  (  ). 

Batchelor  has  shjwn  (see  Hinze,  1959)  that  in  homogeneous  turbulence, 
for  intermediate  times^, 

'hJ  V/j  ,  (  tG) 

i.  e. ,  that 

r  =  /?  e ,/s  ,  ( 


where  $  is  a  constant  and  £  is  the  dissipation  rate  (cm^/scc^). 
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Oamidov  (I960)  suggests  that  ^  is  approximately  equal  to  0.1.  As  will 

-2.  2  3 

be  seen  below,  £  10  /cm  /sec  in  near  surface  waters  of  the  ocean 

_2  -2 

with  moderate  wave  heights.  Assuming  that  €  =10  and  c  =  10"  , 

it  is  found  that  Cv  0.05.  Throughout  the  remainder  of  this  report  it 


will  be  assumed  that 


-  o.  OS’ 


(&) 


therefore 


7(  =  0‘0s  e'h  J 


*h 


($V 


Measurements  of  the  energy  dissipation  rate,  £  ,  as  a  function  of  O 

and  N  thus  yields  information  about  the  corresponding  variation  of  . 

It  is  desired  to  use  equation  (59)  to  calculate  the  diffusion  coefficient 
as  a  function  of  the  energy  dissipation  rate,  £  ,  which,  in  turn,  is  a  func¬ 

tion  of  depth  and  Vaisala  period.  The  validity  of  the  4/3  law  has  been  demon¬ 
strated  experimentally  for  diffusion  in  near  surface  waters,  and  the  constant 
in  equation  (59)  has  been  obtained  from  this  data.  In  order  to  apply  this  rela¬ 
tion  to  intermediate  depths,  it  is  necessary,  at  least,  to  show  that  the  flows 
have  similar  character  and  that  the  application  of  (59)  is  consistent  with  the 
limited  data  available. 


Ozmidov  (1965)  has  suggested  that  a  stable  stratification  can  have  the 
effect  of  suppressing  turbulence  (at  least  the  vertical  motions  of  the  turbulence). 


It  is,  of  course,  well  known  that  in  a  stratified  fluid  with  shear,  there  exists  a 
critical  Richardson  number,  R*  ,  such  that  if 

Ri  ~  ~(f) K*,  (**> 

the  flow  is  unstable  to  small  disturbances  and  presumably  an  unstable  disturb¬ 
ance  will  lead  to  transition  to  turbulence,  (see  Miles  &  Howard,  1964  or,  for  a 
comprehensive  review  of  the  stability  of  stratified  flows,  Drazin  and  Howard,  1966  ). 
Miles  &  Howard  suggest  that  1/4.  If  is  the  flux  Richardson 

number,  see  Lumley  &  Panofsky  (1964),  then  *(  -  (eddy  conductivity/ 

eddy  diffusivity)  1.  Experimental  results  can  be  interpreted  as  giving 
critical  values ,  '/V  a  Rf  £  /. 


Ozmidov  argues  that  if 


e  -  u'/J 


(<u) 


where  it  is  the  turbulent  intensity  and  Jt  is  a  turbulent  scale,  then 

&u/*j  '  uv 
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Lumley  (1964)  has  also  suggested  that,  in  a  stratified  flow,  there 
exists  a  buoyancy  subrange  as  well  as  an  equilibrium  subrange  and  that  the 
wave  number  separating  these  ranges 


-A  =  cj"  *,’/*£-  * 


(it) 


where  C.  is  a  constant  of  order  unity.  If 


V'5/?‘*c% 


(6«y 


(7o; 


The  equilibrium  subrange  of  the  turbulence  exists  for  wavenumbers. 


A.^  A  <  l 


(  /!) 


where  JL  is  the  Kolmogorov  microscale 


A.  -  v- *€  * , 


(7  XJ 


with  Y  the  kinematic  viscosity. 


K  -4.  *-4 


the  equilibiium  subrange  disappears  and  possible. 


there  can  exist  only  a  random  field  of  internal  waves.  For  even  small  scale 
turbulence  to  exist,  there  would  have  to  be  an  increase  in  the  shear  coi-res- 
ponding  to  an  increase  in  £  ,  so  that  would  become  larger  than 

This  general  picture  is  in  agreement  with  the  observations  of  Woods  St  Wiley 
(IS  2)  which  show  that  the  appearance  of  turbulent  patches  is  due  to  short 
internal  waves  in  the  microstructure,  giving  an  increase  in  the  local  shear. 

The  critical  dissipation  rate  for  the  existence  of  an  equilibrium  sub¬ 
range  can  be  estimated  by  equating  *L  and  » 


A/^e;'h  = 


(73) 


(7VJ 


and  for  water  {  I$J  10”^cm^/sec)  and  10  £  10"^sec-^ 


£c  7c  / O  7  cmx / fee7 . 


(75) 
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These  considerations  suggest  an  interesting  experiment.  Produce 
a  turbulent  f'.eld  in  a  stratified  fluid,  perhaps  by  allowing  flow  past  a  grid. 


Initially  £  will  be  large  and  an  equilibrium  range  will  exist.  If  £  is 
large  enough  and  the  experimental  facility  (tunnel)  is  not  very  large,  j( t 

b 

may  be  larger  than  the  diameter  of  the  tunnel  and  a  buoyancy  subrange  may 
not  exist.  Downstream  from  the  grid  the  turbulent  energy  and  also  £  will 
have  decreased;  a  buoyancy  subrange  should  then  exist.  Further  downstream, 
£  will  become  less  than  £  and,  if  this  model  is  correct  the  three- 
dimensional  character  of  the  turbulence  will  entirely  disappear  leaving  a  ran¬ 
dom  field  of  2-dimensionai  turbulence,  possibly  a  random  field  of  internal 


waves . 


It  would  also  be  interesting  to  calculate,  say  from  Ko's  model,  the 
point  at  which  the  energy  dissipation  rate  in  the  wake  drops  below  £^ 

Typical  values  of  ~  C  ~  /  )  are  shown  below.  Near  surface 

O  A 

values  of  £  a^e  of  the  order  of  10  dropping  to  10"*  at  depths. 


Jtc  ('”•) 


^  (cm2/secJ) 

/\/  {se7*"S  «■» 

10-2 
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10’2 

100 

32 

10 

10-3 

3200 

1000 
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It  is  clear  that  a  strong  thermocline  ( 


/J~/o' 


)  and/or  small  energy 


dissipation  rate  will  strongly  suppress  three-dimensional  turbulence.  There¬ 


fore,  the  motion  at  intermediate  depths  is  likely  to  be  two-dimensional,  at 


least  for  scales  larger  than  a  few  meters.  However,  this  is  also  the  case  near 


t.\e  surface  where  the  vertical  extent  of  the  eddies  will  be  of  the  order  of  the 


distance  to  the  surface.  Since  most,  ocean  diffusion  experiments  were  carried 


on  at  depths  of  only  a  few  meters,  the  large  scale  motions  are  likely  to  have 


been  two-dimensional  in  these  experiments  also.  Thus  there  is  at  least  no  in¬ 


consistency  in  the  character  of  the  motions  at  intermediate  depths  and  near  the 


surface. 


4.  2  Experimental  Results 

Woods  &  Fosberry  (1967)  have  made  visual  observations  of  dye 
motions  in  the  thermocline  fine  structure.  They  state  that  ,  the  mini¬ 

mum  eddy  size  observed  varied  from  about  1  mm  to  10  cm,  depending  on  the 
depth  and  stability.  If  ^  is  taken  to  be  the  Kolmogorov  micro-scale  then 

€  =  y3/>?  t  (7t) 

and  corresponding  to  1mm  £  'X  10  cm, 

/o  ,0  c*mX/sec*  £  €  £  /0*a  c^r  */sec\ 

Hale  (1971)  has  observed  visually  the  spread  of  the  dye  above,  in, 
arid  below  a  thermocline  in  fresh  water.  He  states  that  the  turbulent  motion 
appears  to  have  maximum  clergy  near  the  surface  where  it  is  3 -dimensional. 

The  turbulent  intensity  decreases  with  depth;  in  the  thermocline  the  motion 
is  largely  2-dimensional  and  is  perhaps  not  turbulent  at  all.  Below  the  thermo¬ 
cline  the  intensity  increases  again  and  the  motion  becomes  3-dimensional  once 

-I 

more.  Diffusion  coefficients  measured  near  the  surface  were  stated  to  be  un¬ 
related  to  those  at  depth.  I 

Webster  (1969)  has  made  horizontal  current  measurements  at  a  site  ; 

on  the  continental  slope  (39°20'N,  70°  W)  about  175  km  north  of  the  mean  axis  j 

of  the  Gulf  Stream.  The  water  depth  is  2600  m  and  measurements  were  made 

-g 
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from  near  surface  (6  m)  to  a  depth  of  about  2000  m.  He  found  that  some 
of  the  spectra  of  kinetic  energy  density  decreased  as  ^  ^^for  ^  >  0.1 
cycle /hour.  Assuming  that  this  was  an  inertial  subrange,  Webster  calculated 
the  energy  dissipation  rate  as  a  function  of  depth  and  found  that  the  £  thus 
calculated  was  proportional  to  N.  Webster  calculated  values  of  g 

5  x  i0-4cm2/sec3  at  =  10m,  10"4  em^/sec3  at  />  =  100m  and 
10  cm  /sec  at  D  -  1000  m.  A  complete  listing  of  Webster's  calculated 
values  of  £  is  given  in  Table  1. 

Pochapsky  (1972)  measured  the  energy  densities  of  vertical  and 
relative  horizontal  fluid  velocities  with  neutrally  buoyant  floats.  The  measure¬ 
ments,  made  at  a  numher  of  locations,  yielded  spectra  of  the  horizontal  velo- 

»  2 

city  having  a  frequency  variation  between  and  +■  .  Pochapsky 

argued  that  most  of  the  energy  in  this  frequency  range  (0.  1  to  10  eph)  was 
associated  with  internal  waves,  not  turbulence. 

Grant,  Stewart  &  Moillet  (1962)  have  reported  measurements  of 
turbulence  specira  in  near  surface  waters  of  a  tidal  channel.  From  each  of 
their  individual  spectra  they  calculated  energy  dissipation  rates.  These 
ranged  from  0.  02  to  1.0  cm  /sccJ  with  an  average  of  about  0.  34  cm^/sec5. 
These  results  are  listed  in  Table  1. 

Stewart  and  Grant  (1962)  have  reported  measurements  of  at 
various  depths  from  1  to  15  m  and  with  various  wave  heights  in  the  0.  1  to 
0.9m  range.  The  energy  dissipation  rate  showed  a  tendency  to  decrease 


TABLE  I  (cont) 


Data  of  Grant,  Stewart  &  Mpillet  (1962) 

Measured  in  a  tidal  channel  -9  independent  measurements 

£  (cm^/sec^) 

0.610,  1.02,  0.396,  0.121,  0.  235,  0.  147,  0.044,  0.  0187,  0.441 
Data  of  Stewart  and  Grant  (1962) 

Depth  Wave  Height  ^ 


(m) 

(m) 

(cm2/sec3) 

1.0 

0.5 

4.  2  x  10‘2 

1.5 

0.2 

1.  5  x  10'2 

1.5 

0.4 

2.3 

2.0 

0.  1 

5.  2  x  10‘3 

2.  0 

0.3 

2.9  x  10'2 

2.  0 

0.  5 

2.  2 

2.0 

0.9 

4.5 

12.0 

0.  5 

2.  5  x  10'4 

15.0 

0.  1 

1.  1  x  10'3 
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TABLE  I  (cent) 


Data  of  Grant,  Moillet  &  Vogel  (1968) 


Depth 

(m) 

N2 

(sec-2) 

f  3 

(cm2/ sec3) 

1 

esn 

(cm  /sec3) 

15 

1.  56  x  1C"5 

2.  5  x  10"2 

1.00 

2.  5  x  iO"2 

27 

1.56 

5.  2  x  10"3 

1.  00 

5.  2  x  10"3 

43 

1.  56 

3.0 

0.77 

2.3 

58 

5.80  x  10"4 

4.8 

0.  77 

3.7 

73 

2.  05 

1.9 

0.53 

1.0 

89 

1.24 

1.  1 

0.31 

3.4  x  10"4 

90 

1.  24 

1.0 

0.  31 

3.  1 

90 

1.24 

4.8 

0.  31 

1.5 

'wl|lww,a,twlw*Wiinwwiwl^IMPM^*1^^  - ir-TTrrf«rwniirw»rnp^"iiP)w»itirn-  . . . ™ 


with  increasing  depth  and  decreasing  wave  height.  Values  of  £  ranged 
from  4.  5  x  10  ^  to  2. 5  x  10“^  cm^/sec^.  The  values  are  given  in  Table  1. 

Grant,  Moillet  and  Vogel  (1968)  have  measured  £  at  a  number  of 
depths  in  and  above  the  thermocline.  They  found  that  £  decreased  with 

increasing  depth  and  that,  while  the  water  above  the  thermocline  was 
essentially  always  turbulent,  the  turbulence  was  intermittent  in  the  thermo¬ 
cline  and  tended  to  die  out  with  increasing  depth.  Figure  (5)  shows  the  vari¬ 
ation  with  depth  of  the  density  profile,  the  energy  dissipation  rate  (:  , 

the  inter mittency  factor  and  the  "average"  energy  dissipation  rate  fc: 

taken  from  the  paper  of  Grant  et  al  and  Af*  calculated  from 
the  density  profile.  It  can  be  seen  that  £  differs  by  at  least  a  factor  of  5 
across  the  thermocline.  Numerical  values  are  given  in  Table  1. 

Finally  there  is  the  measurement  of  the  diffusion  coefficient  at  300  m 
depth  by  Schu  *rt  (1970).  Schuert  fitted  his  experimental  results  with  a  theo¬ 
retical  model  .  >  *>oint  source  diffusing  according  to  the  4/3  law  and  found 

that  the  measar-  I  value  of  £  was  approximately  1/10  of  the  values  obtained 
in  surface  experiments.  Schuert  found  £  =  6.  6  x  10  cm  *  /sec,  corres¬ 

ponding  to  £  =  2.3  x  10“^  cm^/sec^  (  A  -  0»<?^JLt  a  depth  of  300m  where 
/V«  4  x  10'^sec'^,  (calculated  from  Figure  1  of  Schuert's  paper). 


i 
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It  is  to  be  expected  that  the  energy  dissipation  rates  in  the  tidal 
channel,  as  measured  by  Grant,  Stewart  &  Moillet  (1962)  would  be  con¬ 
siderably  larger  than  those  typical  of  the  open  ocean.  The  current  in  the 
channel  was  several  knots  where  the  velocity  in  the  open  ocean  might  be 
expected  to  be  about  1/4  knot.  For  this  reason  the  tidal  current  measure¬ 
ments  may  be  considered  the  probable  upper  limit  to  £ 

Near  surface  diffusion  experiments,  as  well  as  turbulence  measure¬ 
ments,  are  likely  to  be  made  under  conditions  of  light  winds  and  small  to 
moderate  seas.  The  measurements  of  Stewart  &  Grant  (1962)  are  probably 
representative  of  ’'normal"  conditions.  A  representative  value  of  £ 
near  the  surface  would  appear  to  be  It  cm  /sec  . 

To  summarize,  values  of  £  measured  at  depth  tend  to  lie  in  the  range 

10*  -  10~**  cm  /sec  ,  corresponding,  using  equation  (59),  to  values  of  (5 

•  3  •  3  2/3 

in  the  range  5  x  10"  -  2  x  '10”  cm  '  J/sec.  This  is  in  reasonable  agreement 

3  2/3 

with  Schuert’s  measured  value  of  £  =  6.  6  x  lO"5  cm  / sec  at  a  depth  of  300 
meters.  These  values  of  &  are  about  10  times  smaller  than  that  observed  for 
near  surface  experiments. 


5.  0  Summary  and  Discussion 


I 

I 

1 

T 


The  theoretical  considerations  and  limited  experimental  data  dis¬ 
cussed  in  the  previous  sections  suggest  that: 


1. 


For  a  patch  scale 
dimensional . 


J>Jt 


(equation  67)  the  diffusion  is  two* 


2.  The  vertical  diffusion  coefficient 

-  o,  os  e  J  v/,for  J?  ,  (yj) 

74  *o.ose'/’S//,!0'J>£  (78) 


ft 


3.  The  horizontal  diffusion  coefficient 

X*  0.0s  (7V 

for  Jt  larger  than  the  microscale  and 
smaller  than  the  energy-containing  eddies. 


4.  The  energy  dissipation  rate  £  ,  varies 

with  depth  and  Vaisala  frequency.  Approxi¬ 
mate  values  £  and  the  corresponding 


values  of 

for  different  depths 

are 

€  (cm^/sec^) 

T  (cm2/3/sec) 

Near  surface 

10'2 

10-2 

At  ri.oderate 

** 

1 

0 

1 

CO 

1 

0 

5xl0'3  -  2x  10'3 

depths 


-43- 


i mmtmm 


asusis 


The  calculations  and  results  presented  and  discussed  in  the  previous 
sections  are  based  on  a  simple  model  of  the  diffusion  process  in  the  ocean. 
There  are  a  number  of  important  phenomena  which  have  not  been  included  in 
this  model. 

The  calculations  are  of  the  mean  (ensemble  average)  concentration 
as  a  function  of  space  and  time.  In  any  particular  realization  of  this  field 
there  will  be  statistical  variatiom.  from  the  mean.  The  larger  eddies  will  be 
nearly  frozen  over  periods  of  time  of  the  order  of  ..  For  example, 

turbulent  velocities  will  be  of  the  order  of  lOcm/soc  and  for  eddies  of  the 
scale  of  1/2  km,  the  time  scale  will  be  of  the  order  of  1-1/2  hours.  The 
frozen  or  nearly  frozen  eddies  will  cause  meanders  of  Jie  wake  which  cannot 
be  predicted  from  this  model.  This  model  can  predict  the  variation  of  the 
concentration  off  the  axis  of  the  wake  but  not  the  meanders  of  the  wake  axis. 
See  Schuert  (1970)  for  an  example  of  the  meandering  drift  of  a  single  patch. 

The  model  does  not  include  the  effect  of  gcostrophic  currents  or 
current  shear  on  the  diffusion,  although  the  current  velocities  (a  few  cm/sec) 
are  comparable  to  the  turbulent  velocity.  Nor  have  the  effects  of  inertial 
waves  which  produce  fluid  velocities  comparable  to  that  of  the  turbulent  flow 
field  been  included  in  the  model"  [Pochapsky  (1972)  has  measured  fluid  velo¬ 
cities  due  to  inertial  waves  in  the  range  0.  5  to  10.0  cm/sec|. 

Finally  and  perhaps  most  importantly,  this  model  does  not  include 
the  intermittcncy  of  the  turbulence.  It  is  observed  that  in  a  thermocline 
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turbulence  appears  intermittently,  i.e.,  in  patches.  Woods  &  Fosberry 
(1967)  and  Hale  (1971)  both  comment  on  this  patchiness.  The  measurements 
of  Grant,  Moillet  &  Vogel  (1968)  clearly  show  the  intermittent  character  of 
the  turbulence  with  the  intermittency  factor,  %  £  0.3  at  D  -  90  meters. 

Very  recently.  Woods  &  Wiley  (1972)  have  suggested  a  mechanism 
for  mixing  in  the  interior  of  all  statistically  stable  flows.  This  mechanism 
is  based  on  dye  studies  of  the  microstructure  in  the  thermocline.  Woods  & 
Wiley  call  this  mechanism  billow  turbulence  and  describe  it  as  ", . . .  free 
shear  turbulence  modified  by  a  density  gradient  and  initiated  by  Kelvin  - 
Helmholtz  instability". 

The  Woods  &  Fosberry  mechanism  can  be  described  as  follows: 


The  thermocline  consists  of  thin  sheets,  which  have 
relatively  large  density  changes  across  them  and  more  uniform 
layers  1-2  meters  thick.  There  is  a  velocity  shear  across  the 
layers  due  to  the  geostrophic  current.  Internal  gravity  waves 
propagate  along  the  sheets.  It  is  observed  that,  when  the  sheets 
are  about  10  cm  thick,  short  (  X  ‘>*’75c:n),  steep  internal  waves 
propagating  on  the  sheet  produce  enough  additional  shear  to  make 
the  layer  dynamically  unstable  (  ^  /^**  )•  The  crests  "roll 

up"  and  produce  a  sequence  of  turbulent  billows .  Neighboring 
billows  coalesce,  driven  by  the  shear,  producing  a  turbulent 
patch  of  the  order  of  3  to  5  meters  in  horizontal  extent  and  ap¬ 
proximately  1-2  sheet  thickness  (10  -20cm)  in  vertical  extent. 

The  patch  grows  somewhat  by  entrainment  and  the  turbulence  be¬ 
gins  to  decay.  The  interior  of  the  patch  is  relatively  homogeneous 
but  there  are  large  density  gradients  at  the  top  and  bottom  edges 
of  the  patch.  That  is,  there  are  now  two  very  thin  sheets  at  the 
top  and  bottom  of  the  patch.  These  sheets  spread  by  molecular 
diffusion  until  they  are  about  10  cm  thick  and  then  the  process 
repeats.  This  mechanism  produces  patches  of  fine  scaled,  de¬ 
caying  homogeneous,  isotropic  turbulence.  This  mechanism  is 
assumed  to  be  operative  and  controlling  the  formation  of  turbulence 
throughout  most  of  the  ocean,  i.e.  ,  outside  of  the  near-shore  areas 
and  away  from  the  major  ocean  currents. 
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If  this  mechanism  does  describe  the  generation  and  distribution  of 
turbulence  in  the  ocean,  and  the  evidence  that  it  is  an  accurate  picture  of, 
at  least  part,  of  reality  is  impressive,  then  our  model  of  turbulent  diffusion 
in  the  ocean  must  be  drastically  revised.  The  diffusion  model  would  be  only 
appropriate  as  an  average  over  a  relatively  long  time  scale  (long  compared 
to  the  period  of  patch  gestation) . 

A  more  accurate  model  would  include  large  scale  advection  of  the 
passive  scale  by  internal  and  inertial  waves  ana  ''two-dimensional1'1  turbulence 
(if  it  exists).  At  a  finer  scale,  a  sparse,  random  dist^it>u*ion  of  homogeneous, 
isotropic  turbulent  "pancakes"  would  stir  up  the  scalar  locally.  There  must 
be  a  gap  of  about  one  decade  in  the  energy  spectrum  of  the  fluid  velocities 
since  the  shortest  internal  waves  with  appreciable  energy  (in  this  model)  have 
a  wavenumber 


max 


tjf - x  O.o  /  cm 

75 


,  (* o 


and  the  largest  turbulence  scaler  have  a  wavenumber 


min 


~  O./  CM’'  .  (  fV 

/  o 


A  number  of  approaches  to  this  model  incorporating  large  scale  advection 
and  localized  random  mixing  are  possible.  A  description  of  these  approaches 


will  be  the  subject  of  later  reports, 
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APPENDIX  A 

The  Diffusion  Equation  for  a  Passive  Scalar 

The  discussion  oi  diffusion  given  here  is  based  upon  the  general  dis¬ 
cussion  of  this  phenomena  given  by  Landau  &  Ufshitz  (195?)  but  is  adapted  to 
the  case  of  particular  interest,  namely  diffusion  in  the  ocean. 

Consider  a  mixture  of  two  fluids.  The  equations  of  mass  conservation 
(continuity)  and  momentum  balance  (Navier-Stokes)  remain  unchanged. 


^  "  A  -  °  . 

'  <>x,-  " 


C*-fJ 


where 


^  is  the  total  density, 

(A,  •  is  the  ^  component  of  the  velocitv. 

^  is  the  pressure, 

4A  is  the  1st  shear  coefficient  of  viscosity, 
is  the  2nd  or  bulk  coefficient  of  viscosity, 
is  ^e  ^  ^  space  coordinate,  and 
^  is  the-  time. 
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If  the  relative  concentration  of  the  two  fluids  in  a  given  volume  were 
constant,  the  molon  would  simply  distort  the  volume  and  would  be  absolutely 
reversible.  The  total  change  in  concentration 


OS 

D-t  '  °J 


(*’3) 


with  S  the  concentration  i.  e.  the  ratio  of  the  mass  of  one  component  to  the  total 
mass.  Using  (a-1)  this  is  equivalent  to 

wfrt  +  Txif fsu0  c  °  •  la'VJ 


When  there  occurs  an  interchange  of  material  on  the  molecular  scale,  as 
is  almost  always  the  case,  there  is  an  irreversible  change  in  concentration.  The 
rate  of  change  of  concentration  then  depends  on  the  flux  of  material  at  the  molecu¬ 


lar  level. 


Assuming  that  such  transfer  occurs,  the  transport  equation  for  S  must  be 


modified. 


After  some  algebra,  one  finds 

m/»s  «> 

f(3 r+uJ  <V' 7} 

fKr)^L  .  / 


1  ds  ll 

'  i  r  (*-*) 

\4*/  &  f  Jj  J 
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with 


c/  the  diffusion  coefficient, 

Kt  the  thermal  diffusion  ratio, 

Kjf  Pressure  diffusion  ratio,  and 
j  the  temperature. 

The  second  and  third  terms  on  the  right  hand  side  of  (a- 5)  give  the  change 
in  concentration  due  to  thermal  and  pressure  diffusion.  In  the  cases  of  interest 
there  is  no  measurable  difference  in  diffusivities  due  to  differences  in  physical 
or  chemical  composition.  Therefore 

KT  =  A>  =  o  (a'6j 

to  the  degree  of  approximation  considered  here.  In  this  case, 

-  f  dS  ..  &S  1 


)  t  (i 

, /is ) 

I  ( t 

d  ) 

This  reduces,  in  the  case  of  a  stratified  fluid  (density  grariient  in  the  3  direction  ) 

to  ^  _ 


e>S  ,,.£S  _  ,  O  5  j_J_f  ,)dS 

PT  ■*  /  dYj  '  d  r 

If  the  relative  change  in  desnity  in  the  3  direction  is  small,  and  the  fluid  is  homo¬ 
geneous  in  a  horizontal  plane  (1-2  plane)  then  (a-8)  reduces  to 


ds  j>s_  /  £s_ 

V  ~  d 


L«-9J 
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Now  assume  that  the  fluid  velocity,  ££  *  ,  has  a  mean  (ensemble)  and 


fluctuating  part,  i.  e. 


U:  ~  Uy  +  uj-  t  (O.-/0J 

where,  ^  denoting  an  ensemble  average. 


<  ui>  =  uj 


(a-n) 


Assuming  that 


5 ~  5  +s 


(  a-rtj 


with 


<.s>  -  S  , 


(  a  ~'3J 


equation  (a-9)  becomes,  after  taking  averages  and  assuming  the  statistical  in¬ 


dependence  of  (X/  and  S 

€ 


dS 


...  ds  .  ,2i.\  /  dxs,_ 

7  ^  '  \u/  ■/*  J  ayY}  - 


~'YJ 


Invoking  an  "eddy  diffusivity"  or  mean  gradient  hypothesis 


<.  t  C )s  \  ny '  05 


is 
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equation  (a- 14)  becomes 


5?  +  — 

Sf  1 


•  c)*5“  ^  t  c)S 

d  &px7  +%/  -  ta"CJ 


If  the  fluid  is  stratified,  there  is  no  motion  in  the  3  direction, 
equation  (a- 16)  can  be  written  in  the  form 


c)  5  ,  d5 

*+  +  °i  *% 


)'  e  <>y/  (X/7^1  (Q"V 


wun  / 

*/  ‘J  +  */  ■ 

Assuming  that  the  mean  flow  is  zero, 

u/s  °' 

this  reduces  to  equation  (I). 


(a-  -  /&) 


APPENDIX  B 


Consider  the  equation 


= -J-  fj  Vs  te. 


with  the  initial  jondition 


df 


-), 


(b-0 


SC f,  o)  =  l  t  -'/z  ±f  *  <A  •  Cb-i) 


First,  the  fundamental  solution  will  be  found.  Let 


f  - 


(  b  ~~s) 


Under  this  transformation  equation  (b-i)  becomes 


<^S  /  Z  \  dS 


•(&-)&*(+)  (i-*> 


■7  if/y^ 


Because  (b-4)  differs  only  in  the 


term  from  the  standard 


form  of  the  diffusion  equation,  it  i»  expected  that  its  fundamental  solution  will 
be  similar  to  that  of  the  usual  diffusion  equation. 


Therefore,  let 


S*fy'Vs[  f(ft]  e~<? 


(b-s) 
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with  lo  and  £  constants  and  m  an,  as  yet,  undetermined  function. 
Substituting  this  ansatz  into  equation  (b-1),  it  is  found  that 


J  \  yffSL  n 

In  order  for  this  to  hold  for  «//^  and  ^  , 


(b-Q 


b  =  ¥/?  , 

f  «  yt' . 


(b-7) 

(b-f) 


Thus 


It.  is  clear  that  if 


F* 


(b-/o) 

*h~*o  J  1 


ar.d 


A***-  c  # 


£.  y V*t> — • 


(  b-n) 


showing  that  has  the  character  of  a  & 


function . 


It  is  also  necessary  that 


J **(f,?Vf  =  .»  • 


Equation  (b-12)  is 


(-if)  fe-('**/fl)4-*.  M 


Let 


.  fir-)* 


l  V'SfJ 

and  equation  (b- 1 3)  becomes 


f. 


(I,-/*) 


ilk)  Se-**J*  =  i, 


so  that,  with  (Gradshteyn  &  Ryzhik,  1965,  3.326) 


; 


o 


*-6)r(*l  a-«) 


c  =  i/^nr . 


(b-O) 


With  the  initial  condition  (b-2),  can  be  easily  found  by 


integrating  S  over  the  interval,  thus 


sit*)- 

~'n 


¥ir’AfA  J.,k 


7 


A 


(b-ff) 


tinder  the  transformation 

Z=3(T'<L)'/3Af\  (b-,9) 

it  *  3(f-£/A/zf  >A  /  (b-i  o) 

2z  -  3 (j +i)'h/l  Y,x  .  (&-*') 


Then,  with 


,  "1^4 1,1 1  "  ,fli^,|‘ ' 1 1 1  '  -“i1  T  ,!■  7  f»  vwirniHW*!  yen^jiiff  1,'Wwwn  jf. ft  ftfi  % *•  “ 


- . - r~Trpf>',*rff'ff>  n  ■~rip(TVTTtt’i''’t‘IJ<((;]^.i.fu*~ «w«^H 


S(f,  V  -  if  e*i(*r)  -  C'-ffa) 


fP 


[z,e-^-?xe-*']J  (b.zi) 


The  concentration  at  the  center  of  the  patch  is  given  by 


S(o,p-rf[jiy/,]  -[J^Je-hA'^ Ui)  j 


As  for  the  scale  or  size  of  the  patch.,  it  can  be  defined  as  in  the  case  of  the 
constant  eddy  diffusivity, 

S C'k't  y)  ~  &*  S{ 0,  y)  ;  Ib-W) 

with  Q  c  a*<f.  Equation  (b-24)  can  be  solved  numerically  for  ^  as  a 
function  of  ^  ,  using  (b-22)  and  (b-23). 


Defining  the  variance  as  usual  by 


rz=  J  f'Sfcn)  <ff. 


(b-zrj 


Defining  the  variance  as  usual  by 


r’* 


(  b-ZC) 


it  can  be  seen,  substituting  from  equation  (b-18)  that. 


r 


With  the  substitution 


and  an  interchange  in  the  order  of  integrations 


J  '{ Ijr) 


where. 


m)  ?« is  *  (sbh  v'x*fc 


ib-”j 


l y" e'*  ■  a'n) 


Carrying  out  the  integrations 


( b-yfj 
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